Vortex Lattices in Rotating Atomic Bose Gases with Dipolar Interactions 
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We show that dipolar interactions have dramatic effects on the groundstates of rotating atomic 
Bose gases in the weak interaction limit. With increasing dipolar interaction (relative to the net 
contact interaction), the mean-field, or high filling fraction, groundstate undergoes a series of tran- 
sitions between vortex lattices of different symmetries: triangular, square, "stripe", and "bubble" 
phases. We also study the effects of dipolar interactions on the quantum fluids at low filling frac- 
tions. We show that the incompressible Laughlin state at filling fraction v — 1/2 is replaced by 
compressible stripe and bubble phases. 
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Ultra-cold atomic Bose gases have emerged as remark- 
able systems with which to study the unusual response of 
Bose-condensed systems to rotation. Experiments have 
allowed the imaging of lattices of quantised vortices 0, 
and the study of their collective dynamics^- New vor- 
tex lattice structures in two-component condensates have 
been observed Q, and a novel regime of vortex density 
in which the vortex cores overlap [j] has been accessed 
Theory shows that, at very high vortex density 0], atomic 
Bose gases should undergo a transition into novel uncon- 
densed phases closely related to the incompressible liq- 
uids of the fractional quantum Hall effect 0, 0, 0, 0- 

Although Bose condensation can occur for a non- 
interacting Bose gas, the formation of arrays of quan- 
tised vortices under rotation relies on non-vanishing (re- 
pulsive) interparticle interactions. In typical atomic 
Bose condensates, the interactions are so short-ranged 
that they can be viewed as local (contact) interac- 
tions. However, significant additional non-local inter- 
actions can arise if the atoms have intrinsic or induced 
electric or magnetic dipole moments [!£ [I0j. The recent 
achievement 1 1 lj of the Bose condensation of chromium 
(which has a large permanent magnetic dipole moment) 
has opened the door to the experimental study of dipolar- 
interacting Bose gases. It is important to ask what are 
the effects of non-local interactions on the properties of 
the vortex lattices. 

In this paper, we study the effects of dipolar interac- 
tions on the groundstate of a weakly interacting atomic 
Bose gas under rotation. We show that the additional 
non-local interaction leads to dramatic changes in the 
nature of the groundstate. Within mean-field theory, the 
triangular vortex lattice predicted for contact interac- 
tions is replaced by vortex lattices of different symme- 
tries: square, "stripe" and "bubble" phases. Further- 
more, we study the properties at very high vortex den- 
sity where the groundstates for contact interactions are 
incompressible quantum fluids. We show that with in- 
creasing dipolar interactions the Laughlin state, which 



is the groundstate for pure contact interactions at filling 
fraction v = 1/20, is replaced by compressible states 
that are well described by stripe and bubble phases. 

We consider a system of bosonic atoms of mass M con- 
fined to a harmonic trap with cylindrical symmetry about 
the z-axis. We denote the natural frequencies of the trap 
by u>\i and luj_ in the axial and transverse directions, and 
the associated trap lengths by ay j_ = d%] (Muj\\_±). The 
particles are taken to interact through both contact inter- 
actions and additional dipolar interactions. We consider 
the atomic dipole moments to be aligned with the z-axis, 
as in the experiments reported in Ref. 0] , such that the 
net two-body interaction is 

where a s is the s-wave scattering length arising from con- 
tact interactions an d Cd is a measure of the strength 
of the dipolar interactions. 

We study the regime of weak interactions, when the 
mean interaction energy per particle is small compared 
to the trap energies, fkj±». The single particle states 
are then restricted to the groundstate of the trap in the 
z-direction (2D) and to the lowest Landau level (LLL) 
of the x — y motion^], 0]. While non-rotating gases are 
typically far from the weak interaction limit, this limit 
can be approached at high angular momentum owing to 
the reduction in particle density through the centrifugal 
spreading of the cloud[j|. In this limit the interactions 
are fully specified by the Haldane pseudopotentialsQ, 
V m : the interaction energy of a pair of particles with 
relative angular momentum m. For bosons, symmetry 
of the wavefunction means that only even m contribute. 
The precise form of the dipolar pseudopotentials depends 
in a non-trivial way on the trap asymmetry a^/aj_. For 
simplicity of presentation, we study the limit a\\/a± — > 
0, in which the dipole forces fall off most quickly with 




FIG. 1: Contour plots of the particle densities for condensed states on a torus with Ny = 16 vortices. The light (dark) shading 
indicates high (low) particle density (on arbitrary scales), (a) Triangular vortex lattice (a/b = v3/2); (b) Square vortex lattice 
(a/b = 1.0); (c) "Stripe" phase (shown for a = 0.528, a/b = 0.608); (d) q = 4 "bubble" phase {a/b = a/3/2). 



increasing distance. To leading order in a\\/a±, we find 
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Vo represents the net Zoca/ interaction, which has contri- 
butions from both the contact and dipolar interactions; 
V m >o represent non-local interactions which arise only 
from the dipolar interaction. We consider the relative 
sizes of the non-local and local interactions to be vari- 
able, either by tuning the dipolar interaction, Cd, or by 
tuning a s close to a Feshbach resonance^) (which can 
allow a s to become negative). We quantify their relative 
size by the ratio 
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We are interested in bulk properties, so perform calcula- 
tions on a periodic rectangular geometry (a torus) which 
accurately describes the centre of an atomic gas where 
the particle density is approximately uniform. In the 
weak interaction limit the number density of vortices is 
ny = l/(7raj_)j so a torus with sides a and b contains 
Ny = a6/(7Tflj_) vortices. 

First, we describe the mean-field groundstates. Within 
Gross-Pitaevskii (GP) mean-field theory the groundstate 
is assumed to be fully condensed, with condensate wave- 
function, i/>(r). I n the weak-interaction limit, ip(r) is 
found by minimising the mean interaction energy for 
fixed average particle density, with the wavefunction con- 
strained to states in the 2D LLL^j|. The problem 
is mathematically equivalent to the Ginzburg-Landau 
model for a type-II superconductor close to H c2 01 gen- 
eralised to a non-local interaction. A phenomenological 
model of this kind has been discussed in the context of 
superconductivity ^3|> but we are not aware of general 
solutions for the vortex lattice groundstates. We have 
found the mean-field groundstates for the interactions 



by numerical minimisation on a torus with up to 
Ny = 24 vortices (i.e. exploring periodic states with 
up to 24 vortices in the unit cell). As a function of a, 
the groundstate undergoes a series of transitions between 
states of different translational symmetries. Representa- 
tive images of the particle distributions in these states 
are shown in Fig 2] 18] . 

The groundstates we find are: a triangular lattice of 
single vortices (0 < a < 0.20); a square lattice of sin- 
gle vortices (0.20 < a < 0.24); a stripe phase (0.24 < 
o < 0.60) 19]. The stripe phase consists of broad lines of 
high particle density separated by rows of closely-spaced 
vortices. The vortices are ordered along the rows, so the 
state is a "stripe crystal" with crystalline order in both 
directions. For a > 0.60 the states consist of clusters 
of high particle density arranged in a triangular lattice. 
Owing to the similarity to crystalline states of electrons 
in high Landau levels |20j. we refer to these as "bubble" 
states. We find a sequence of bubble states, which we 
label by the number q of vortices associated to each bub- 
ble. The bubble states we find as groundstates are q = 4 
(0.60 < a < 0.91), q = 5 (0.91 < a < 1.4), q = 6 
(1.4 < a < 2.0), q = 7 (2.0 < a < 2.7), q = 8 (2.7 < a). 
The particle distributions of the bubble states with q > 4 
resemble Fig^d), with additional vortices confined to 
the honeycomb network separating the bubbles of par- 
ticles. We have not looked for states with q > 8, but 
expect that states of arbitrarily large q will appear. 

We believe that it should be possible experimentally 
to access these new vortex lattice groundstates. To make 
contact with experimental parameters, we now examine 
the case of a spherically symmetric trap (an = aj_). In 
this case the triangular lattice is replaced by new ground- 
states if a s is tuned[li to values a s < — 0.l3CdM/h . 
Note that, despite this negative value of a s the contribu- 
tion of the dipolar interaction makes Vq positive, which, 
in the weak interaction limit, is sufficient to ensure sta- 
bility to collapse. Under these circumstances, for a mean 
particle density 713^, the interaction energy per particle 
is of order Cdn^d, so the weak interaction limit requires 



(a) 



0.5 

a 



(b) 



1X1 



13.4' 



(c) 



3 

\K\ 



FIG. 2: Results of exact diagonalisation studies at v — 1/2. (a) Overlap of the exact groundstate with the Laughlin state 
(N = 8, Nv = 16, torus aspect ratio a/b = 1). (b) Energy spectrum for a = 0.528 (N = 9, Nv = 18, a/b = 0.82). The 
filled circles are for momenta K — (0,0) and (0,3), showing evidence for a stripe state, (c) Energy spectrum at a = 0.758 
(N = 8, Nv = 16, a/b = V3/2). The filled circles are at K = (0, 0), (4, 0), (0, 4), (4, 4), (2, 0), (1, 2), (3, 2), (2, 4), which are the 
reciprocal lattice vectors of the q — 4 bubble state. The shaded symbols are at K = (0, 2), (2, 1), (2, 3), (4, 2), corresponding to 
the reciprocal lattice vectors of a rotated crystal. 



n>3d ~ hoj±JCd- Choosing the value of Cd appropriate for 
chromium | ll|. which has a magnetic dipole moment of 6 
Bohr magnetons, and taking loj_ = 2ir x lOOrad s~ , the 
weak interaction limit should be a good approximation 
for densities less than about 10 14 cm~ 3 . 

We now turn to discuss the groundstates beyond the 
mean-field approximation. In RefofiJ it was shown that 
the parameter controlling the validity of mean-field the- 
ory for a rotating atomic Bose gas is the filling fraction, 
v = n/ny, where n and ny are the number densities 
(per unit area) of particles and vortices. For contact 
interactions, the triangular vortex lattice predicted by 
mean-field theory was shown to be destroyed by quan- 
tum fluctuations for v < vf 1 , with i/* n ~ 60, and re- 
placed by new groundstates which include incompress- 
ible liquids closely related to fractional quantum Hall 
states H 0,013. 

We have studied the effects of dipolar interactions on 
these strongly-correlated groundstates using exact di- 
agonalisation studies. We now return to the case of 
Oy/oj. — * and focus our attention on the filling frac- 
tion v = 1/2, for which the exact groundstate for con- 
tact interactions (a = 0) is the v = 1/2 bosonic Laugh- 
lin state 4]. This is the dominant incompressible state 
of rotating bosons with contact interactions. We find 
that the exact groundstate in the presence of dipolar 
interactions is well described by the Laughlin state up 
to a ~ 0.5. At this point there is an abrupt transi- 
tion in the groundstate and its overlap with the Laughlin 
state falls to a very small value [FigEta)]- Our studies 
show that the states that replace the Laughlin state at 
a > 0.5 are compressible phases which are well described 
by the mean-field groundstates discussed above. Here we 
present evidence of a stripe phase at cv = 0.528 and of 
the q — 4 bubble phase at a = 0.758. (For larger a we 
find evidence of bubble phases with larger q.) To iden- 



tify these broken-symmetry states we make extensive use 
of the classification of energy eigenstates by a conserved 
moment um |2 1 \ . We express the momentum as a dimen- 
sionless vector K = [K X ,K V ) using units of 2nh/a and 
2i:h/b for the x and y components, and report only pos- 
itive K x ,K y [states at (±K x ,±K y ) are degenerate by 
symmetry]. We identify the broken translational sym- 
metry of the groundstates by making use of the fact that 
this leads to the appearance of quasi-degeneracies in the 
spectrum at wavevectors equal to the reciprocal lattice 
vectors of the crystalline order [22^. 

Evidence of stripe ordering at a = 0.528 is presented 
in Figj2£b). This shows the excitation spectrum on a 
torus with an aspect ratio chosen to be consistent with 
the mean-field groundstate. The quasi-degeneracy of the 
groundstate at K — (0, 0) with a state at (0, 3) indicates 
a strong tendency to translational symmetry breaking 
in the stripe pattern found in mean-field theory: three 
stripes lying parallel to the short axis of the torus. How- 
ever, unlike the mean-field state, we find no evidence of 
crystalline order parallel to the stripes: the groundstate 
appears to be a "smectic" [ijjj in which translational sym- 
metry is broken in only one direction. 

Evidence for the formation of the q = 4 bubble phase 
is shown in Fig. Etc). At filling fraction v — 1/2, each 
bubble contains vq — 2 particles, so this state can also be 
described as a triangular lattice of pairs of particles. The 
low energy states shown in Fig. Etc) as filled symbols cor- 
respond to wavevectors at the reciprocal lattice vectors 
of the q = 4 bubble state [shown in FigHJd)]. Athough 
these states do not seem to be cleanly separated from the 
rest of the spectrum, the other low energy states (shaded 
symbols) can be understood as arising from a nearby (in 
energy) crystalline configuration in which the crystal is 
rotated through 90°. 

The appearance of stripe and bubble states at v = 
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1/2 indicates that these states are much more stable 
to quantum fluctuations than is the triangular vortex 
lattice (which is unstable for v < v* n ~ 6.6|). This 
enhanced stability can be understood within a simple 
Lindcmann analysis, in which one asserts that quan- 
tum melting occurs when the quantum fluctuations of 
a lattice site exceed a multiple Ci of the lattice spac- 
ing. Treating the fluctuations of each vortex indepen- 
dently, one expects the triangular vortex lattice to melt 
for v < ^* ri = \/3/(27r4)ja For the square lat- 

tice, we find v s ^ = 1/(ttc|) = (2/^)^, close to that 
for the triangular lattice. The enhanced stabilities of the 
stripe and bubble phases arise from the existence of larger 
lengthscales and larger numbers of particles per unit cell 
(as compared to the triangular or square vortex lattices). 
For the q bubble phase, the bubbles form a triangular lat- 
tice with lattice constant ^Jq^a± . Applying the Lindc- 
mann analysis to the quantum fluctuations of the centre- 
of-mass of a bubble, we find = ^4-4 = z/'"/? 2 - Even 

for the smallest (q = 4) bubble state, this critical filling 
fraction is very much smaller than that for the triangular 
lattice. For the stripe phases, there are two lengthscales: 
the inter- vortex separation in the directions parallel, Ru , 
and perpendicular, R±, to the stripes. One therefore ex- 
pects two transitions: when fluctuations of the vortices 
along the stripes exceed clR\\ there is a loss of order in 
that direction, leading to a smectic phase; when fluctua- 
tions perpendicular to the stripes exceed c^R± the stripe 
ordering will finally be lost. Assuming that, for this fi- 
nal transition, of order R±/R» vortices must fluctuate 
together, we find that loss of stripe order should occur 

at = ^ (^) 2 = (2/V3) (R l{ /R ± )\f. Over 

the range a = 0.24 — 0.60 for which the stripe is the 
mean-field groundstate the ratio R±/R\\ = 1 — 2.59, so 
the critical filling fraction for the stripe can be as small 
as about (l/6)z/* n . While it would be desirable to have 
a more complete Lindemann analysis in which the col- 
lective modes|25J of these new lattices are quantised, we 
believe that the simple analyses presented here capture 
the essential physics of the relative stability of the states 
to quantum fluctuations. 

At filling fractions above v = 1/2 quantum fluctua- 
tions of the stripes and bubbles are strongly suppressed. 
For the cases in Fig. Gib) and (c), increasing the num- 
ber of particles to N = 12 (so that the respective filling 
fractions are v = 2/3 and 3/4) leads to much improved 
groundstate quasi-degeneracies. Quantum fluctuations 
are enhanced for v < 1/2. For small non-zero a, we find 
incompressible liquids at filling fractions, v = p/ (3f>± 1), 
expected for composite fcrmions 8] formed from bosons 
bound to three vortices. The state at v = 1/4 is well de- 
scribed by the v — 1/4 Laughlin state. For a > 1.7 this 
state is replaced by bubble phases, consistent with the 
expectation from the Lindemann analysis that for large 



q these can have critical filling fractions less than 1/4. 
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